Universal decay of classical Loschmidt echo of neutrally stable but mixing dynamics 
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We provide analytical and numerical evidence that classical mixing systems which lack exponential 
sensitivity on initial conditions, exhibit universal decay of Loschmidt echo which turns out to be a 
function of a single scaled time variable <5 2//5 i, where 8 is the strength of perturbation. The role of 
dynamical instability and entropy production is discussed. 
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Fidelity, or Loschmidt echo, is defined as the over- 
lap of two time evolving states which, starting from the 
same initial condition, evolve under two slightly differ- 
ent Hamiltonians. It is therefore an important quantity 
which measures the stability of the motion under systems 
perturbations. The recent interest in the behaviour of fi- 
delity [JSSSSESEIlllilinihas been largely 
motivated by a possible use in quantifying stability of 
quantum computation Il2| |. 

It has been shown [lfj that for classical chaotic, ex- 
ponentially unstable systems, the decay rate of fidelity 
is perturbation independent and, asymptotically, fidelity 
decays as correlation functions. On the other hand, for 
quantum systems, fidelity decay obeys different regimes 
depending on perturbation strength. In this relation, 
particularly intriguing is the recently discovered case of 
mixing dynamics with vanishing Lyapounov exponent 
[T3I a prominent example of which are billiards inside 
polygons [ijj . In several respects the classical dynamics 
of such systems is reminiscent of quantum dynamics of 
generic chaotic systems which, apart from an initial time, 
exponentially short in h, are linearly stable. As a conse- 
quence statistical relaxation in quantum mechanics takes 
place in absence of exponential instability. Certainly, the 
dramatic difference in the dynamical stability properties 
of different systems must be reflected in a different qual- 
itative behaviour of physical quantities such as fidelity 
which is the object of the present paper. 

In the following, under the assumptions of linear sepa- 
ration of trajectories and dynamical mixing, which may 
be produced by some discontinuity in the flow, we de- 
rive a universal scaling law of classical fidelity decay. We 
conjecture that this surprising fidelity decay may be as- 
sociated to a peculiar power-logarithmic entropy produc- 
tion in such systems. We consider here a specific exam- 
ple, i.e. the triangle map z n+ \ = T(z n ) |14| on a torus 

z = (x,y) e [-1,1) x [-1,1) 



y n+1 = y n + asgna;„ + (3 (mod 2), 
x n+ i = x n + y n+1 (mod 2), 

where sgnx = ±1 is the sign of x and a, (3 are two 
rameters. Previous investigations have shown that 



(1) 



(see also 0] for some rigorous results on {Q): for ra- 
tional values of a, (3 the system is pseudo-integrable, as 
the dynamics is confined on invariant curves. If a = 
and (3 is irrational, the dynamics is (uniquely) ergodic, 
but not mixing, while for incommensurate irrational val- 
ues of a, {3 the dynamics is ergodic and mixing with dy- 
namical correlation functions decaying as i~ 3 / 2 . It can 
be argued that the triangle map possesses the essential 
features of bounce maps of polygonal billiards and Id 
hardpoint gases 0, namely parabolic stability in 
combination with decaying dynamical correlations, and 
as such represents a paradigmatic model for a larger class 
of systems. 

The classical fidelity Fs(n) can be written as an over- 
lap of two phase space densities propagated by the orig- 
inal map T and the perturbed map T$ = To g$ where 
gs(z) — z + 5a(z) is some near-identity area-preserving 
map parametrized by a vector field a(z): 



F 5 (n) = 



Jd 2 zp(T(-Kz))p(T^ n, (z)) 
J d 2 zp 2 (z) 



(2) 



We can make our discussion even more general by tak- 
ing the perturbation explicitly time- dependent. Let the 
perturbed map Tg iTl explicitly depend on iteration time, 
namely we consider the following class of perturbed tri- 
angle maps, z n+ i = Ts,n(z n ) 

Vn+1 = Vn + asgnx n + [3 + Sf n (x n ) (mod 2), 
x n+ i = x n +y n+ i (mod 2). (3) 

We will assume that the force function f n has vanish- 
ing time-average for almost any initial condition. Let as 
further assume that the initial density p(z) is a charac- 
teristic function over some set A of typical diameter oj 
with 5 <C oj <C 1. Then a pair of orbits z n and z n start- 
ing from the same point zq = zq in A, contribute to @ 
until they hit the opposite sides of the discontinuity, at 
x = 0, 1 (mod 2). The fidelity at time n is then simply 
the probability that the pair of orbits does not hit the 
cut up to n— th iterate. Assuming ergodicity of the map 
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FIG. 1: Fidelity decay. To compute fidelity we divide the 
phase space in 100 x 100 cells. We then take 10 points in 
one cell and evolve these points with the map Q up to time 
n. Then we compute the reverse evolution with the map 
© and compute the fraction of points which fall again in 
the initial cell, after time 2n. The result is then averaged 
over initial distributions in 49 different randomly chosen cells, 
(a) log(l — Fs(n)) versus log((5 2 ' 5 n) magnifying behaviour for 
short times. The broken curves refer to S = 10" 5 , . . . , 10" 10 
(top-down). The full line is the theoretical expression ©. (b) 
log[(F<s(ra) — u)/{l — u)] versus 5 2 ^ 5 n magnifying asymptotic 
(long-time) behavior. Here u = 10~ 4 is the relative area of 
the initial set A giving the asymptotic value of fidelity. The 
meaning of broken curves is the same as in (a). The full 
curve gives the numerical solution of the random Gaussian 
model 11121 . while the dashed line has slope —1 to indicate 
asymptotic exponential decay. 



Il4l we write 



F B (n) = ( A (1 - |Aav 



(4) 



\n' = l 



where Ax n = x n — x n , and (A n ) = 
J dzp{z)A{T^ n \z)) / J dzp(z). In order to derive 
the fidelity decay for the triangle map we have to 



compute the average growth rate of the orbits distance 
perpendicular to the cut. This is achieved by writing out 
an explicit linearized map for the orbits displacement 

Az n — z n z n 



Az„ + i = 



1 1 
1 



Az n 



5fn{x n ) 



(5) 



This system of linear difference equations can be solved 
explicitly, say for Ax n : 



Ax n 



n'=0 



(n - n')f n >(x„>) 



(6) 



with the initial condition Azq = 0. Assuming that f n 
are pseudo-random variables with quickly decaying corre- 
lation function C(n) = lim m _,.oo (fm(%m)fn+m(Zn+m)), 
we can employ a version of the central limit theorem to 
show that Ax n should have Gaussian distribution for suf- 
ficiently large n. To this end, let us first notice that the 
second moment 



n— 1 n — 1 
n'—0 n"=Q 



can be related, as n — > oo, to the integrated correlation 
function. Since, for large n: (fn'fn") — C(n' — n"), we 
obtain by means of a straightforward calculation 



<(Az„) 2 ) 
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C{m). (8) 



Now, as long as fidelity remains close to 1 , we can expand 
Q to first order Fg(n) = 1 - YZ>=1 (\ Ax n>\) , where 
the average (|Ax n |) = (3tt)|(5||77.| 3 / 2 can be com- 

puted using a Gaussian distribution of Ax n with variance 
((Ax„) 2 ) given in JSJ. This yields 



F s (n) 



757T 



\S\\n 



5/2 



(9) 



This expression is valid until F$(n) remains close to 1, 
that is up to time |n| < n* = x / 5 |<5| 2/5 

In fig. la we show the behavior of 1 — Fg(n) for short 
times n < n* and compare with the theoretical formula 
@ with a = 3.29 ± 0.01 as computed from numerical 
simulation of correlation function C(n). As for pertur- 
bation we choose a simple shift in the parameter a, so the 
force reads f(x) — sgnx and is, in this case, not explicitly 
time dependent. Yet it is pseudorandom and one can see 
that, as S decreases, the numerical curves approach the 
theoretical expression ©. 

Notice that according to eq. JSJ, the average distance 
between two orbits increases as cx n 3 / 2 . On the other 
hand, the distance between two initially close orbits of 
the same map increases only linearly with time. This is 
nicely confirmed by the numerical simulations of fig. 2. 

For larger times, n > n* , higher order terms in the ex- 
pansion of 10} contribute, so temporal correlations among 
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FIG. 2: (a) Average distance d n — \Az n \ versus time n for 
two nearby initial orbits of the unperturbed map The ini- 
tial distance is do = 10~ 9 , the average is taken over 2.5 x 10 8 
different initial conditions. The dotted line has slope 1. (b) 
Average distance versus time for two orbits starting from the 
same initial condition and evolving under the unperturbed 
and perturbed maps, Q and ©, respectively. The val- 
ues of perturbation 5, for the curves (from top down), are: 



10" 



10" 



The data are averaged over 



10 J different initial conditions. Dotted lines have slope 1.5. 




FIG. 3: Spreading of phase space points of the echo- 
dynamics. We consider 5000 initial points in the central cell 
of a 3 x 3 grid (a). We then evolve these points up to time 
n, and then reverse the motion with the perturbed dynamics, 
with 8 — 10~ 6 , up to the echo time 2n. The density of points, 
for n = 200, 400, 600 is shown in (b),(c),(d), respectively. 



Ax n become important. We are here unable to de- 
rive exact theoretical predictions for the fidelity decay 
in this regime. However, numerical results in fig. lb 
show that, for large times, fidelity decays exponentially 
Fs(n) = cxp(— 7|n|) with exponent 7 = C|<5| 2 / 5 . We also 
checked that the transition time between the two regimes 
of decay, scales as S~ 2 ^ 5 . In conclusion, extensive and 
accurate numerical results provide clear evidence that fi- 
delity depends on the single scaling variable r = 5 2 / b n. 

In the following, we show that this scaling behavior 
can be derived analytically for sufficiently small 6. The 
only assumption is correlation decay with a finite char- 
acteristic time-scale n m ix, i-e. (fnfn') practically vanish 
for \n — n'\ > w m j x . Let us divide the time-span n into 
v := n/m blocks of m steps each, such that n m i x <C m <C 
n„ and make a scaling argument. The local variation of 
Ax n , namely Ax n+ i-Ax n = <5En'=o/"' ~ <5\/nismuch 
smaller than the mean value (|Ax„|) ~ Sn 3 / 2 . Thus we 
approximate the product within each block labelled 



by i = 1, . 
Therefore 



v as (1 |Ax (t _ 1)m |) r 



1 



\Ax 



Fs{n)^l[(l-m\Ax {i _ 1)m \). 



(i-l)ml 



(10) 



Next we define the normalized block-averaged forces 




-l)m+k 



(11) 



k=0 



which are normalized, and uncorrelated, = <5 tjU 

since m ^> n m ; x . Using Eq.© we can write Axu-i) m ~ 

^E^=l(i-M)"*Efc^o'/(M-l)"»+fc = <fr™ 3/2 0- 1/2 E^=l(i- 

/i)£p. If, in additon to the rescaled time v — n/m, we 



define a rescaled perturbation e 
write Eq. (fTU|) as 



&a x l 2 m?l 2 then - 



*«(")= (n i 1 - 



(12) 



The derived relation Fs(n) = & Sa i/2 m 5/2(n/m) does not 
depend on m (for large enough m) , and therefore fidelity 
should be a function of the scaling variable r = |<5| 2 / 5 n 
only. 

Notice that due to the central limit theorem, since 
m 3> Umixj £fi can be simply treated as uncorrelated, 
normalized, Gaussian stochastic variables. We have ac- 
tually computed the universal function </>(e 2 / 5 ^) = $ e (^) 
by means of Monte-carlo integration, and checked that it 
is practically insensitive to e, for e < 10~ 4 . As it is seen 
in fig. lb, the numerical data for the triangle map agree 
with the theoretical expression (|12fl . namely 0((5 2 / 5 cr 1 / 5 n) 
which is plotted as a full curve. 

The two regimes of fidelity decay described above are 
illustrated in fig. 3 by the image at the echo time of an 
initial uniform phase space distribution over some set A. 
Notice that the linear-response regime (0 is valid until 
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FIG. 4: The time evolution of the coarse-grained entropy 
for the triangle map, computed by taking 5 x 10 7 points ini- 
tially distributed randomly in one cell of a TV x TV phase-space 
grid with TV = 700. We plot ln(S(n) - 5(1)) versus In Inn. 
The straight line has slope 3/2. Obviously, the entropy will 
eventually saturate at S(oo) = In TV 2 . In the inset we show 
the entropy computed for the uniquely-ergodic and nonmixing 
dynamics with a = 0. The dotted line has slope 1. 

the shape of the intial set is approximately restored at the 
echo time. For larger times, the fidelity decay becomes 
exponential. 

Finally we would like to stress that this behavior of 
the triangle map differs from the typical behavior which 
has been found for chaotic or for integrable systems. In 
particular, contrary to the case of exponentially unstable 
systems, in this case the rate of fidelity decay depends 
on the perturbation strength. This feature is shared by 
quantum systems in which exponential instability is ab- 
sent as well. One may wonder if this behavior is re- 
flected also in some other, perhaps even more funda- 
mental dynamical property of the map. In order to ex- 



plore this question, we have computed the entropy pro- 
duction for the map £[J. As the extensive computation 
of Kolmogorov-Sinai dynamical entropies seemed too ex- 
pensive for reaching any conclusive results, we have de- 
cided to compute the dynamical evolution of the coarse 
grained statistical entropy S n — —J2jPn^ n Pn^- To 
this end we divide the phase space in TV x TV equal cells, 
and consider an initial ensemble of points uniformly dis- 
tributed over one cell. The probability Pn is defined 
as the fraction of orbits which, after n time steps, fall 
in the cell of label j. For a chaotic system with dy- 
namical entropy h, one expects S n = hn + const |l7j . 
whereas for ergodic-only (non-mixing) dynamics one ex- 
pects S n ~ Inn, for sufficiently large TV. Our numerical 
results for the triangle map (fig. 4) show instead that 
S n — Si = | lnn| A with the exponent A = 3/2. Further- 
more, as shown in the inset of fig. 4, for the triangle map 
with a — numerical results give, quite accurately, 
S(n) = Inn (with no prefactor or additional constant). 

In conclusion, we have discussed the parametric sta- 
bility, as characterized by classical fidelity or Loschmidt 
echo, of an important class of dynamical systems where 
neutral stability is coexisting with dynamical mixing. As 
a paradigmatic example of this class of systems we have 
considered the triangle map. By means of analytic cal- 
culations and numerical simulations we have derived two 
universal regimes of fidelity decay, both being character- 
ized by a universal scaled time variable |5| 2 / 5 i. This inter- 
esting dynamical behavior is supported also by a power- 
logarithmic behavior of the coarse-grained entropy. 
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